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ABSTRACT 

In a previous work [I], we proposed an optimal multi-path selec- 
tion method for multiple description (MD) encoded video stream- 
ing. 1 To do so, we first modelled multi-path streaming and then 
developed an expression, i.e., an objective (cost) function, that es- 
timated average streaming distortion in terms of network statistics, 
media characteristics and application requirements. Naturally, the 
ultimate goal was to find the set of paths that minimized this cost 
function. However, finding such sets of paths turned out to be in- 
tractable in large topologies. Hence, in this paper, we provide a 
fast heuristics-based solution by exploiting the infrastructure fea- 
tures of the Internet. The simulations run over various random 
Internet topologies show that the proposed heuristic is able to find 
a good solution in a much shorter time than the brute-force ap- 
proach. Particularly, this heuristic is best suited to such interactive 
multimedia applications as video-conferencing and VoIP, where 
multi-path computation is a time-critical process. In addition, it 
is also suitable for the clients whose processing power capabilities 
are limited. 

1. INTRODUCTION 

Multiple Description Coding (MDC) is a source coding that gener- 
ates several encoded sub-bitstreams called descriptions, which are 
transmitted between the server and client via partially link-disjoint 
paths. These descriptions can still be decoded independently. At 
the client, if all descriptions are received error-free, then a maxi- 
mal quality signal reconstruction is possible. However, if at least 
one description is available at the client, then the client can still 
reconstruct a signal of acceptable quality. MDC provides this ro- 
bustness at the expense of some reduction in the compression effi- 
ciency. The reduced correlation between packet loss probabilities 
on partially link-disjoint paths makes MDC a good choice for loss 
resiliency. Moreover, in contrast to ARQ and FEC, MDC not only 
provides error resiliency features, but it also avoids excessive de- 
lays; thus, MDC has been found to be an effective solution for 
interactive multimedia applications in lossy networks [2, 3]. 

Several studies have proposed multi-path transport (path di- 
versity) to improve the reliability of the streamed video [3, 4, 5]. 
Also there are some studies to model path diversity based on the 
loss characteristics of the paths in point-to-point and multipoint- 
to-point networks [6, 7]. However, these studies do not propose a 
formal path selection method that considers all aspects of network 
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1 It should be mentioned that the proposed multi-path selection method 
is optimal within the context of the network models and objective functions 
that are considered in [1]. 



statistics, media characteristics and application requirements. To 
address this problem, we modelled multi-path streaming and pre- 
sented a formal way to evaluate the streaming distortion for each 
set of paths in a given network [1]. Consequently, the optimal se- 
lection was a straightforward evaluation of distortion for each pos- 
sible set and choosing the one with the minimum distortion. How- 
ever, the complexity of this approach increases exponentially in 
the total number nodes. Therefore, optimal computation may not 
be feasible if the client does not have enough processing power, 
memory space and time to find the solution or if the client resides 
in a network where the conditions rapidly change. To overcome 
these problems, we investigate fast heuristics in this paper. For the 
sake of continuity, we first give a brief overview of the multi-path 
modelling. Interested readers are referred to [1] for further details. 

2. MULTI-PATH STREAMING MODELLING 

Let's consider two paths, V\ and V2, between the nodes Ns and 
Nc, which are referred as a path pair. Our ultimate goal is to find 
the pair that minimizes the average distortion. For this purpose, 
we need to estimate the end-user quality in terms of network pa- 
rameters. In case of two descriptions, we can write the average 
distortion at the client as 

Po,0 

D = Prob {Both received on time} x £>o,o (1 ) 

^0,1 

' * H * 

+ Prob {1 st received on time & 2 lost or delayed} x D 0 , 1 
+ Prob {1 st lost or delayed k 2 nd received on time} x Z)i,o 

^1,0 

H- Prob {Both lost or delayed} x£>i,i . 

Here, D m ,n and P m , n (m, n G {0, 1}) denote a specific distortion 
and its occurring probability, respectively as defined in Table 1 . 
In order to capture the characteristics of multi-path streaming, we 
seek a simple yet representative model. Fig. 1 depicts our end-to- 
end multi-path model. In this figure, we have two paths between 
the server and client that split at node Nz> We combine the joint 
links between the server and Nz to form a sub-network denoted by 
Networkjomt. On Networkj 0 i„f, the two correlated descriptions 
are supposed to be transmitted back-to-back and hence, they are in 
the same burst period. On the other hand, we consider the rest of 
the paths as disjoint despite a possible merge at a later node. That 
is, the two correlated descriptions do not fall in the same burst pe- 
riod implying that their loss events beyond Nz will be independent 
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Fig. 1. End-to-end multi-path model 



of each other. This is a valid approximation as long as the disjoint 
links have a delay difference larger than the average burst period 
between them [1]. This is typically true because of the unevenly 
experienced delays and jitters on disjoint links. Thus, the links be- 
yond Nz are combined to form the respective Network^ S j oint . 
Note that node Nz might be the server, which implies that the 
paths are totally link-disjoint, i.e., Networkjoint is an empty set. 
The pertaining parameters to our model are tabulated in Table 1 . 
We start our derivations with the distortion. 

In compressed-video applications, the distribution of Discrete 
Cosine Transform (DCT) coefficients can be modelled as Lapla- 
cian. For a Laplacian source with variance a 2 , the distortion as a 
function of the aggregated bandwidth of the path pair (B) can be 
approximated by [8] 



D(B)~ a 2t T 2 ^ 



(2) 



where W x H (pixels/frame) and F (frames/second) represent the 
frame resolution and frame rate, respectively, c is a known con- 
stant depending on the chroma subsampling format. Given this 
relation, we can compute the distortion terms in (1) by plugging 
the corresponding B values into (2). First, for Do,o the aggregated 
bandwidth is typically the sum of the individual bandwidths, i.e., 
B = B\ + Bi provided that the joint links have available band- 
widths higher than this sum. However, if this is not the case, B is 
then limited by the bottleneck joint-link capacity. Second, to com- 
pute £>o,i and Di,o, we take B — B\ and B = B2, respectively. 
Finally, D\,\ is computed by letting B = 0. 

The next step in the computation of (1) is to find Po,o — 
each of which is comprised of two components. The first compo- 
nent, arrival probability, computes the probability that the descrip- 
tion is not lost during the transmission and the second component, 
on-time delivery probability, computes the probability that the de- 
scription arrival occurs before a pre-specified deadline. 

We adopted a well-known two-state Markovian Gilbert-Elliott 
(GE) model to describe the temporal behavior of packet losses on 
a link. We give the transition probabilities, a and /?, for the link 
l u>v by 



Pti,v X fiu t V 

1 — Pu,v 



and 



fiu,v — j 



(3) 



We employ a joint GE model for Networkj oiTl t and individual 
Bernoulli models for Networki isjoint based on our model. To 
find the joint GE parameters, namely a joint and Pjoinu we first 
approximate the joint packet loss rate and average burst length for 
Networkjoint as follows [1]: 



Vi u>u € Network Joint 



Lj = — x 
Pj 



(4) 



(5) 



l>u,v 

b u ,v 
Vi 

Dm.n 



Bm,n 

Bi 

TTA 

a 

Llu,V 

Pj 

Jj 

A 
Ji 



Table 1. Notation and model parameters. 



Node k 

Link between the nodes N u and N v 

Bandwidth on link l UyV 

Packet loss rate on link l u , v 

Minimum (constant portion) delay on link Z«, v 

Jitter on link l U}V 

Pathi 

Distortion, m and n denote the arrivals of 1 st and 2 nd 
descriptions, respectively. 

0 — > arrival on time. 1 — » loss or excessively delayed. 

Probability corresponding to D m ,n 

End-to-end bandwidth on path Vi 

Steady-state prob. for state A (arrival), tta = 

Steady-state prob. for state L (loss), ttl = 

Transition prob. from A to L 

Transition prob. from L to A 

Average burst length on link 

Packet loss rate on Networkj 10int 

Average burst length on Network joint 

Packet loss rate on Network^ oint 

Jitter on Network joint. 

Jj = Dvi u>v 6Network ioint 

Jitter on Network^ oint . 

3d. = SvZ UfV 6Networkj, iajoiTit 3n,v 

End-to-end jitter on path Vj. Jj = jj + j d 



Then, we compute a joint and 0 joint by plugging pj and Lj 
into (3). On the other hand, since they are modelled by Bernoulli, 
the packet loss rates in Network^ oint are given by 

Pd-i- n (i-p« (V ). (6) 

Vi«, v 6NetworkJ itfioint 

The second step in computing the probability P m , n is to derive 
the on-time delivery probabilities. Usually, different paths have 
different tolerances for jitter. We define t^ aa: as the maximum jit- 
ter that can be tolerated for path Vi. LetT* = £ V j u v€V . t u , v rep- 
resent the end-to-end minimum delay on path Vi. Then, we have 
tUax = T dead u ne - Ti where T dead u ne is the tolerable delay for 
the target application. Now, the on-time delivery probability can 
also be given as the probability that J» is smaller than t^ax. 

Finally, we have all the information required to write success 
probabilities (Po,o — Pi,i) explicitly. The corresponding equa- 
tions are given in (7) - (10) where we denote the joint GE pa- 
rameters by a and (3 instead of a joint and Pjoinu respectively. 

Po.o =((tta(1 - a) 2 + 7r LJ 9(l - a)) x (1 - p l d ) x (1 - p 2 )) (7) 
x Prob{jj + j d < t l max and jj + j d < t 2 max ) 



V/„, v e Networkj- oint 



Po,i =((ir A (l-ot)a + -x L 0a)x (1 - pj)) 
xProbUj 

+ (0m(1 - a) 2 + * L 0(1 - a)) x (1 - pi) x p 2 ) 
x Prob(jj + j l d < t l nax ) 

+ ((7^(1 - a) 2 + * L P(1 - a)) x (1 - pi) x (1 - p 2 )) 
x Prob(jj + j d < tl ax and jj + j 2 > &„) 



(8) 



Pi.o =((7r A a/3 + 7r t (l-/3)jS)x(l-p=)) (9) 
xProbUj+jl <t 2 max ) 

-f ((;r A (l - a) 2 + 7^0(1 - q)) x x (1 - p5)) 
xProbUi+fi <t 2 max ) 

+ ((^(1 - a) 2 + tt l /?(1 - a)) x (1 - pi) x (1 - p 2 )) 
x ProbUi + ii > tl ax and + jj < t 2 nax ) 

=l-Pb,o-Po,i -Pi.o (10) 

3. SOLUTION APPROACHES 

Once we have defined the cost function, we seek the optimal path 
pair such that the expression in (I) is minimized. Given that (1) is 
highly nonlinear, solving this optimization problem may become 
computationally intractable in large topologies. In this section, we 
study the brute-force complete enumeration and a heuristic that 
selects a good pair of diverse paths. This heuristic is based on 
first identifying a small subset of good paths, and then choosing 
the pair that minimizes (1) among this subset. In developing our 
heuristic, we exploit the routing hierarchy in the Internet. See [9] 
for further details on this subject. 

3.1. Brute-Force Approach 

Consider an Internet topology with M transit domains (TD) each 
of which contains N nodes. An illustrative topology with four 
TDs is shown in Fig. 2. Assume that multi-path routing capability 
is enabled only within TDs, but not stub domains (SD). In this 
topology, the brute-force (BF) approach for finding the optimal 
pair enumerates all possible paths between the server and client. 
Then, (1) is computed for each pair and the pair that gives the 
minimum average distortion is selected. 

The complexity analysis of this approach can be summarized 
as follows: All paths in the given topology can be generated in 
0(2 MN ) time and (I) is evaluated for all possible pairs in 0(MNx 
2 2 N ) time. Finally, identification of the optimal pair requires 
0(log2 2MN ) time. Hence, the BF approach has a complexity of 

0(2 MN )^(MNx2 2MN )-K)(MN)==0(MNx2 2MN ). (11) 

In (11), we observe that the BF approach complexity increases 
exponentially with both M and N. Below, we develop a more 
efficient heuristic that is referred as rapid path generation (RPG). 

3.2. Rapid Path Generation 

An important characteristic of the Internet is that a path connecting 
two nodes in different SDs, e.g., the ones containing the server 
and client, will never pass through any other SDs [9]. This feature 
implies that multi-path routing will mainly take place within or 
between the TDs. Moreover, between two TDs, usually there is 
only one link through which all traffic from the source domain is 
carried to the destination domain. This natural decomposition of 
the Internet leads us to consider the TDs individually in the path 
generation process. 

We start our algorithm by finding good paths in each TD. Dur- 
ing this process, we also exploit the fact that among the links in a 
TD, only a few of them may be heavily congested at a time while 
the rest is lightly congested or not congested at all. Eliminating 
these congested links clearly reduces the topology size. Moreover, 
after this reduction the remaining links will have similar (almost 




Fig. 2. An illustrative Internet topology with four TDs. 



zero) packet loss rates. On the other hand, the end-to-end band- 
width is often limited by the SD or inter-domain links. We suppose 
that TD links can support the bandwidth desired by the clients. 
Hence, the remaining TD links' bandwidths have negligible effect 
in selecting the good paths. 

Now, we are left only with the delay and jitter parameters in 
the path generation process. To find a subset of good paths in each 
TD, we solve a A:-shortest path algorithm where the cost on each 
link is the sum of its delay and average jitter. The most straight- 
forward fc-shortest path algorithm for M TDs requires 0(kMN 2 ) 
time. After selecting k good paths in each TD, we combine them to 
create 0(k M ) end-to-end paths. Then, out of these paths we form 
0(k 2M ) pairs. At this point, we evaluate (1) for each pair. This 
process takes 0(MNx k 2M ) time. Finally, selecting the pair with 
minimum average distortion can be accomplished in 0(logk 2M ) 
time. Hence, the RPG approach has an overall complexity of 

0(kMN 2 )+0{MN x k 2M )+0(M) = 0(MN x k 2M ). (12) 

Generally, we have M < N, which is the case in small-to- 
moderate sized networks. Hence, although the RPG approach runs 
in exponential time in M, it still provides an important improve- 
ment over the BF complete enumeration, which requires exponen- 
tial time in the total number of nodes (M x N). In addition, with 
this approach rather than trying to seek the solution in the whole 
network, we handle reduced-sized problems in each TD separately. 
Hence, the demand for the memory space in the RPG implementa- 
tion is substantially lower than that in the BF implementation. On 
the other hand, k is an important parameter in (12). We will see its 
impact on the performance of RPG approach in the next section. 

3.3. Performance Comparisons 

In order to evaluate the performance of RPG, we used random In- 
ternet topologies generated by GT-ITM [9]. We varied k, number 
of TDs and SDs to see the effects of variations in network con- 
nectivity. The link parameters were assigned random values as de- 
scribed in [1]. For each random topology, we first found the path 
pairs corresponding to the BF and RPG solutions. Then, by using 
each pair we streamed a standard test sequence from the server to 
the client with a delay tolerance of 200 ms. 2 We repeated this pro- 
cess 1 00 times to obtain reliable peak signal-to-noise ratio (PSNR) 
values. The results are tabulated in Tables 2-4 along with the val- 
ues of k, M and N. The PSNR differences between the BF and 



2 We used Table Tennis sequence. This sequence comprised of 150 
frames each of which was 352 x 240 pixels. The frame rate was 30 f/s. 



RPG approaches are also given in the last column of each table. 
Note that each row represents a different random topology. Hence, 
the results in a row of a table should only be compared with the 
results in the same row of other tables. 

Tab le 2. Comparison of the BF and RPG approaches when k = 1. 



k 


M 


N 


BF 


RPG 


Difference 




2 


5 


33.39 dB 


31.82 dB 


1.57 dB 




2 


8 


33.45 dB 


31.56 dB 


1.89 dB 




3 


5 


33.36 dB 


31.63 dB 


1.73 dB 




3 


8 


33.38 dB 


31.44 dB 


1.94 dB 




4 


8 


33.30 dB 


31.39 dB 


1.91 dB 



Tab le 3. Comparison of the BF and RPG approaches when k = 2. 



k 


M 


N 


BF 


RPG 


Difference 


2 


2 


5 


33.39 dB 


32.93 dB 


0.46 dB 


2 


2 


8 


33.45 dB 


32.85 dB 


0.60 dB 


2 


3 


5 


33.36 dB 


32.84 dB 


0.52 dB 


2 


3 


8 


33.38 dB 


32.66 dB 


0.72 dB 


2 


4 


8 


33.30 dB 


32.43 dB 


0.87 dB 



Table 4. Comparison of the BF and RPG approaches when k — 3. 



k 


M 


N 


BF 


RPG 


Difference 


3 


2 


5 


33.39 dB 


33.32 dB 


0.07 dB 


3 


2 


8 


33.45 dB 


33.30 dB 


0.15 dB 


3 


3 


5 


33.36 dB 


33.24 dB 


0.12 dB 


3 


3 


8 


33.38 dB 


33.15 dB 


0.23 dB 


3 


4 


8 


33.30 dB 


33.06 dB 


0.24 dB 



The results show that as we identify more good paths in each 
TD, the performance of RPG approach becomes closer to the op- 
timal in all topologies. Although there is about 1.80 dB difference 
in the average qualities when k = 1, this difference reduces to 
about 0.16 dB when k = 3. This is also shown in Fig. 3 and 4 
where individual frame qualities are plotted. In Fig. 4 the RPG 
plot tracks the BF plot more closely on the average than it does in 
Fig. 3. Even when k = 3 (and M = 4, N = 8), the number of 
operations required by the RPG approach is 500K times smaller 
than the one required by the BF approach, which proves a quite 
large savings in terms of processing power. 

Another interesting point is the choice of fc. In addition to the 
reported results, we also tested the same topologies with k = 4. 
However, the improvement over the case with k = 3 was not that 
significant. Hence, we conclude that keeping k — 3 is a good 
choice. It is worthy to note that k = 1 does not necessarily imply 
that the 1 st and 2 nd shortest paths are selected between the server 
and client. Although only the shortest path is computed in each 
TD for k = 1, we still evaluate (1) for all possible end-to-end 
pairs and choose the pair accordingly. 

4. CONCLUSION 

In this paper, we presented a fast heuristics-based solution for the 
optimal multi-path selection problem introduced in [1]. This ap- 
proach performs within 0.2 dB of the optimal solution in small-to- 
moderate sized networks without incurring any run-time or mem- 
ory space problems. Particularly, this fast heuristic is best suited 
to such interactive multimedia applications as video-conferencing 
and VoIP, where multi-path computation is a time-critical process. 
In addition, it is also suitable for the clients whose processing 
power capabilities are limited. 
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Fig. 3. The PSNR comparisons of the BF and RPG approaches 
when k = 1, M = 4 and N = 8. 
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Fig. 4. The PSNR comparisons of the BF and RPG approaches 
when k = 3, M = 4 and N = 8. 
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